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Abstract. We present a three-dimensional model of rotating convection combined with a sim- 
plified corona in spherical coordinates. The motions in the convection zone generate a large-scale 
magnetic field that gets sporadically ejected into the outer layers above. Our model corona is 
approximately isothermal, but it includes density stratification due to gravity. 
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1. Introduction 

The Sun sheds plasma into the heHosphere via coronal mass ejections (CMEs). There 
has been significant progress in the study of CMEs in recent years. In addition to im- 
proved observations from spacecrafts like SDO or STEREO, there have also been major 
advances in the field of numerical modeling of CME events. One of the main motiva- 
tions for understanding the generation and dynamics of CMEs is to have more reliable 
predictions for space weather. CMEs can have strong impacts on Earth and can affect 
microelectronics aboard spacecrafts. However, an important side effect of CMEs is that 
the Sun sheds magnetic helicity from the convection zone which may prevent the solar 
dynamo from being quenched at high magnetic Reynolds numbers (Blackman & Bran- 
denburg, 2003). 

In most of the CME models, the motion in the photosphere as well as the initial and 
footpoint magnetic fields are prescribed or taken from two-dimensional observations at 
the solar surface. Such fields represent an incomplete approximation to the full three- 
dimensional field generated by the dynamo. An alternative way of modeling CMEs would 
be to perform a 3-D convection simulation to generate the configuration of magnetic field 
and the photospheric motions self-consistently. However, the convection zone and the 
solar corona have very different timescales. In solar convection the dominant timescale 
varies from minutes to days. While this is short compared with the dynamo cycle, time 
scales in the solar corona can be even shorter because the Alfven speed is large. 

In earlier work (Warnecke & Brandenburg 2010, Warnecke et al. 2011) we have estab- 
lished a two-layer model with a unified treatment of the convection zone and the solar 
corona in a single three-dimensional domain. In those models, the generation of magnetic 
fields was driven by turbulence from random forcing with helical transverse waves mim- 
icking the effects of convection and rotation in a simplified way. This two-layer model 
was able to produce recurrent plasmoid ejections which are similar to observed erup- 
tive features on the Sun. In this work we develop this approach further and apply a 
self-consistent convection model instead of random forcing. This model automatically 
includes differential rotation as a result of the interaction of rotation and convection. 
Here we present some preliminary results from such a study. We find the formation of a 
large-scale magnetic field, which eventually gets ejected as a magnetic structure. It seems 
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to us that this mechanism could play an important role for the formation of CMEs and 
flare events on the Sun and other stars. 



2. The model 

As in Warnecke & Brandenburg (2010) and Warnecke et al. (2011) a two- layer model 
is used. Our convection zone is similar to the one in Kapyla et al. (2010, 2011). The 
domain is a segment of the Sun and is described in spherical polar coordinates (r, 6>,^). 
We mimic the convection zone starting at radius r = 0.7 R and the solar corona until 
r = 1.5 i?, where R denotes the solar radius, used from here on as our unit length. In 
the latitudinal direction, our domain extends from colatitude = 15° to 165° and in 
the azimuthal direction from (/) = 0° to 90°. We solve the equations of compressible 
magnetohydrodynamics. 
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T^ = -V- KVT + 21.5^ + ^ J2 - reooi, (2.4) 
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where the magnetic field is given by B = V x A and thus obeys V-S = at all times. The 
vacuum permeability is given by po , whereas magnetic diffusivity and kinematic viscosity 
are given by r] and respectively. Sij = ^{Ui-j + Uj-i) — l^ij V • U is the traceless rate- 
of-strain tensor, and semicolons denote covariant differentiation, Q = Qo{cosO^ — sin^, 0) 
is the rotation vector, K is the radiative heat conductivity and g = —GMr/r^ is the 
gravitational acceleration. The fiuid obeys the ideal gas law with 7 = 5/3. We consider 
a setup in which the stratification is convectively unstable below r = Rsun^ whereas the 
region above is stably stratified and isothermal due to a cooling term Fcooi in the entropy 
equation. The Fcooi term is r dependent and causes a smooth transition to the isothermal 
layer representing the corona. 

The simulation domain is periodic in the azimuthal direction. For the velocity we use 
stress- free conditions at all other boundaries. For the magnetic field we adopt radial field 
conditions on the r = 1.5 boundary and perfect conductor conditions on the r = 0.7 
and both latitudinal boundaries. Time is expressed in units of r = (lirms^f)"^, which is 
the eddy turnover time in the convection zone. We use the Pencil CODEyj, which uses 
sixth-order centered finite differences in space and a third-order Runge-Kutta scheme in 
time; see Mitra et al. (2009) for the extension to spherical coordinates. 



3. Results 

In this work we focus on a run which has fiuid Reynolds number Re = 3, magnetic 
Reynolds number Rcm = 32 and Coriolis number Co = 7. We define the Reynolds num- 
ber as Re = Urms/^h for the magnetic one Rcm = Urms/v^^U respectively and the Coriolis 
number as Co = 2fl/urmsh' After around 100 turnover times, the onset of large-scale 
dynamo action due to the convective motions is observed. The magnetic field reacts back 
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Figure 1. Variation of B(p and Br in the convection zone at r = 0.89R (left panel) and 
r = 0.79R (right). Dark blue shades represent negative and light yellow positive values. The 
dotted horizontal lines show the location of the equator at = 7r/2. The magnetic field is 
normalized by the equipartition value. 
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Figure 2. Time series of coronal ejections in spherical coordinates. The normalized current 
helicity, fxoRJ - B/ {B^)t, is shown in a color-scale representation from different times; dark 
blue represents negative and light yellow positive values. The dashed horizontal lines show the 
location of the surface at r = i?. 



on the fluid motions and causes saturation after around 200 turnover times. The satu- 
ration is combined with an oscillation of the magnetic field strength in the convection 
zone. The field reaches its maximum strength of about 60% of the equipartition field 
strength, (5eq = (moP'^^)^^^, which is comparable with the values obtained in the forced 
turbulence counterparts both in Cartesian and spherical coordinates (Warnecke & Bran- 
denburg 2010, Warnecke et al. 2011). The magnetic field in rotating convection seems to 
show certain migration properties. In Figure [TJ we show the azimuthal (^0) and radial 
(Br) magnetic fields versus time (t/r) and latitude (90° — 6) for two different heights. 
The magnetic field emerges through the surface and is ejected as isolated structures. 
The dynamical evolution is clearly seen in the sequence of images of Figure [21 where the 
normalized current helicity {jnoRJ • B/ {B'^)t) is shown. If one focuses on the region near 
the equator {0 = 7r/2), a small yellow (i.e. positive) feature with a blue (negative) arch 
emerges through the surface to the outer atmosphere, where it leaves the domain through 
the outer boundary. This ejection does not occur as a single event — it rather shows recur- 
rent behavior. We do not, however, find a clear periodicity in the ejection recurrence, like 
in earlier work. Even though the ejected structures are much smaller than in Warnecke 
et al. (2011), their shape is similar. However, the detection of an ejection with the aid 
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Figure 3. Dependence of the dimensionless ratio /jqRJ • B/ {B'^)t on time t/r and radius r in 
terms of the solar radius. The top panels show a narrow band in 6 in the northern hemisphere 
and the bottom ones in the southern hemisphere. We have also averaged in latitude from 4.1° 
to 19.5° (left panel) and 32.5° to 45.5° (right). Dark blue shades represent negative and light 
yellow positive values. The dotted horizontal lines show the location of the surface at r = i?. 

of the current helicity is much more difficult in convection-driven simulations than in 
forced turbulence. In Figure [21 one can see large structures diffusing through the surface 
into the upper atmosphere at higher latitudes. These structures disturb the emergence 
of ejections and hamper their detection. These larger diffusive structures are also visible 
in Figure [3l where the normalized current density is averaged over two narrow latitude 
bands on each hemisphere. The formation of these diffuse structures in the corona seem 
to get suppressed when the stratification of the system is increased. Nevertheless, ejec- 
tions are still visible, for example around t/r = 1900, t/r = 2200 (see Figure Ej) and 
t/r = 2400. 

The hemispheric rule, as we have found in earlier work, is not completely valid in the 
current work. The current helicity in Figure [3] shows a different behavior in the northern 
and southern hemispheres, but one cannot tell clearly the leading sign. This has to do 
with the much lower values of relative kinetic helicity h^e\{r^ t) = {uj • u) / uj^msUTms in the 
convection simulations. Values of up to h^Q\ = ^0.4 are reached at some radii in the two 
hemispheres. In the forced turbulence simulations of our earlier work, we studied purely 
helical systems with nearly /irei = Tl- 

In summary, we are able to advance our two-layer model approach by including self- 
consistent convection, which generates the magnetic field and eventually drives ejections. 
Not surprisingly, the ejections occur non-periodically and as smaller structures than in 
earlier work, but this behavior seems to be similar to the Sun. Furthermore, detailed 
investigations covering a wider range of magnetic and kinetic Reynolds number as well 
as rotation rates show promising results. 
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